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L  INTRODUCTION 


The  response  of  a  broad  range  of  materials  (ceramics  including  optical  windows,  geophysical 
solids,  polymers,  composites,  etc.)  to  intense  impulsive  loading  is  central  to  many  Air  Force  pro¬ 
grams  and  applications.  With  rapid  advances  in  materials  processing  and  fabrication  technology 
for  novel  materials,  and  the  wide  spectrum  of  applications  envisioned  for  these  materials,  it  is  diffi¬ 
cult  to  test  materials  under  all  conditions  of  interest.  Hence,  a  fundamental  understanding  of 
mechanisms  governing  material  response  under  impulsive  loading  is  needed.  Understanding  and 
predicting  the  highly  nonlinear  material  response  for  a  diverse  range  of  materials  to  extreme  load¬ 
ing  conditions  is  a  challenging  problem.  Knowledge  gained  from  quasi-static  deformation  studies  is 
of  limited  value  because  the  underlying  deformation  mechanisms  at  high  strain-rates  can  be  quite 
different.  In  the  case  of  energy  deposition  from  an  intense  laser  pulse,  there  is  the  additional  need 
to  understand  the  fundamental  processes  controlling  the  rapid  energy  deposition  in  materials  prior 
to  the  so-called  damage  threshold. 

The  overall  objective  of  this  research  effort  was  to  understand  the  nonlinear  material  response 
due  to  very  rapid  energy  deposition  in  selected  solids.  Two  different,  though  related  aspects  of  the 
overall  problem  were  addressed  in  the  present  work.  In  the  first,  the  response  of  brittle  materials 
to  plane  shock  waves  (produced  by  projectile  impact)  was  examined  using  novel  methods.  The 
thrust  of  this  effort  was  to  understand  the  mechanisms  for  material  deformation  at  very  high 
strain-rates  and  stresses.  In  the  second  aspect  of  the  research  program,  the  detailed  microscopic 
mechanisms  that  are  responsible  for  the  absorption  of  laser  energy  in  transparent  dielectrics  and 
subsequent  mechanical  and  thermodynamic  response  of  the  material  were  investigated. 

An  extended  summary  of  the  shock  wave  and  the  laser  deposition  efforts  is  presented  in  Sec¬ 
tions  II  and  III,  respectively.  The  emphasis  in  these  sections  is  on  problem  definition  and  summar¬ 
izing  the  research  accomplishments.  Detailed  accounts  of  these  efforts  may  be  seen  in  journal  arti¬ 
cles  and/or  dissertations  that  are  listed  in  Appendix  A.  These  articles  and  dissertations  have  been 
submitted  to  AFOSR.  A  related  goal  of  the  present  work  was  the  training  of  graduate  students 
and  interaction  with  national  laboratories.  Information  regarding  these  objectives  is  listed  in 
Appendices  B  and  C,  respectively. 
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II.  SUMMARY  OF  THE  SHOCK  WAVE  EFFORT 


Plane  shock  wave  experiments  provide  a  unique  method  for  understanding  material  response 
at  very  high  stresses  and  strain-rates.  The  ability  to  make  time-resolved  measurements,  with 
nanosecond  resolution,  under  well  defined  strain  conditions  is  particularly  attractive  for  studying 
time-dependent  processes  responsible  for  material  deformation.  Despite  the  extensive  literature  on 
the  subject,  there  remain  two  main  difficulties  in  understanding  micro-mechanical  mechanisms 
under  shock  loading:  (i)  the  inability  to  experimentally  determine  the  complete  stress  state  in  the 
material  and  (ii)  the  lack  of  direct  experimental  information  at  the  microscopic  level.  The  present 
work  addresses  these  long  standing  problems  using  recent  experimental  developments. 

The  primary  goals  of  our  work  were  to  examine  the  response  of  selected  materials  using  com¬ 
bined  compression  and  shear  loading,  and  to  explore  the  development  of  time-resolved  optical  spec¬ 
troscopy  for  studying  high  strain-rate  deformation  of  solids.  The  following  four  tasks  comprised  the 
shock  effort:  (i)  Shear  wave  measurements  in  shocked  fused  silica  to  understand  structural  changes 
accompanying  nonlinear  elastic  deformation  in  this  material,  (ii)  Use  of  combined  compression  and 
shear  loading  in  shocked  calcium  carbonate  (Carrara  marble)  to  understand  structural  phase  transi¬ 
tions  and  to  examine  the  interrelationship  between  material  strength  and  phase  transitions,  (iii)  To 
explore  the  development  of  time-resolved  spectroscopic  methods  to  understand  shock  induced 
deformation  in  single  crystals  at  the  microscopic  level,  (iv)  Exploratory  effort  to  examine  shock 
induced  phase  changes  in  calcite  single  crystals.  The  single  crystal  calcite  effort  was  motivated  pri¬ 
marily  by  the  comprehensive  effort  undertaken  in  Task  (ii).  Further  discussion  regarding  these 
tasks  is  summarized  below. 

A.  Nonlinear  elastic  response  of  shocked  fused  silica 

Fused  silica  displays  anomalous  compression  under  both  hydrostatic1,2  and  shock  loading3 
conditions.  The  purpose  of  the  present  work  was  to  use  shear  wave  measurements  to  obtain  a 
more  complete  description  of  the  nonlinear  elastic  response  of  fused  silica  in  the  shocked  state  than 
is  currently  available. 

Using  in-material,  particle  velocity  gauges  in  inclined,  parallel  plate  impact  experiments,4  5 
large  amplitude  compression  and  shear  wave  profiles  were  measured  in  fused  silica  samples  (Dynasil 
1000).  The  longitudinal  stress-strain  relation  calculated  from  the  compression  wave  profiles  was 
found  to  be  identical  over  the  stress  range  (to  60  kbar)  of  the  present  work  to  that  reported  for 
G.E.  Type  15i  fused  silica3 

Shear  wave  velocities,  in  the  shocked  state,  were  found  to  decrease  linearW  with  compression 
up  to  a  compression  of  7.6  percent.  At  higher  compressions,  the  shear  wave  velocities  increased 
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with  compression.  These  data  in  conjunction  with  the  measured  longitudinal  stress-strain  relation 
were  used  to  obtain  the  longitudinal,  shear  and  bulk  moduli  in  the  compressed  state.  An  interest¬ 
ing  finding  w as  that  the  bulk  and  shear  moduli  display  a  minima  at  different  density  compressions. 
In  addition  to  these  moduli,  the  Poisson’s  ratio  for  shocked  fused  silica  was  determined  as  a  func¬ 
tion  of  density  compression.  It  displays  a  complex,  nonlinear  behavior  with  values  ranging  between 
0.17  and  0.06  over  the  12  percent  density  compression  range  examined.  This  finding  contradicts 
the  assumption  of  constant  Poisson’s  ration  for  fused  silica.6  and  has  interesting  implications  for 
stress  deviators  in  shocked  fused  silica. 

The  mean  stress  versus  density  curve  for  shocked  fused  silica  was  calculated  by  integrating 
the  bulk  modulus  versus  compression  curve.  The  mean  stress-density  compression  curve  under 
shock  wave,  uniaxial  strain  loading  differs  appreciably  from  the  measured  hydrostat  and  provides 
direct  evidence  for  shear  enhanced  compaction  in  fused  silica.  This  phenomenon  is  completely 
reversible  to  60  kbar.  The  observed  shear  enhanced  compaction  appears  to  be  compatible  with  a 
model  proposed  by  Hemley  et  al.7  for  the  structural  changes  in  hydrostatically  compressed  fused 
silica.  Definitive  support  of  this  conjecture  will  require  Raman  measurements  under  shock  loading. 

Using  the  mean  stress-density  and  longitudinal  stress-density  relation  obtained  in  the  present 
work,  we  obtain  very  large  values  of  stress  deviators  in  shocked  fused  silica.  The  large  stress  devia- 
tor  values  may  play  a  significant  role  in  the  phase  changes  in  fused  silica  observed  at  higher  pres¬ 
sures.8 

The  present  data  were  used  to  calculate  second  and  third  order  elastic  constants  for  fused  sil¬ 
ica.  Two  third  order  constants  (C1U  and  C112)  were  obtained  from  the  present  data.  Two 
independent  but  consistent  methodologies  were  used  in  analyzing  our  data.  Both  gave  approxi¬ 
mately  the  same  results. 

B.  Shock  induced  phase  transformation  in  polycrystalline  calcite 

Structural  phase  transformations  under  shock  loading  are  of  considerable  interest  for  develop¬ 
ment  of  constitutive  relations  at  high  strain  rates  and  for  understanding  the  response  of  solids 
under  nonhydrostatic  stresses.  Over  the  stress  range  where  both  phase  changes  and  mechanical 
yielding  can  occur,  it  is  difficult  to  separate  these  two  phenomena  from  longitudinal  wave  profile 
measurements.  Because  crystal  structures  cannot  be  directly  measured  under  shock  loading,  it  is 
difficult  to  identify  the  phase  transition  mechanisms.  Comparisons  with  static,  high-pressure  data 
are  needed  to  interpret  the  shock  data.  These  comparisons  are  difficult,  for  materials  with  signifi¬ 
cant  strength,  due  to  the  difficulty  in  separating  effects  due  to  stress  deviators  and  loading  rates. 
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In  the  present  work,  the  response  of  polycrystalline  calcium  carbonate  was  examined  to  assess 
the  usefulness  of  measuring  both  shear  and  compression  waves  for  better  understanding  shock- 
induced  phase  transformations  in  solids.  In-material  particle  velocity  histories  were  measured  in 
Carrara  marble  subjected  to  combined  compression  and  shear  loading.  These  wave  measurements 
presented  the  possibility  of  determining  additional  information  on  the  shear  response  of  the 
shocked  material,  and  also  the  complete  stress  state  produced  by  uniaxial  strain  loading.  Both  pos¬ 
sibilities  could  contribute  to  improved  understanding  of  shock-induced  phase  transformations  by 
permitting  more  direct  comparisons  with  static,  high-prc^ure  data,  and  providing  greater  informa¬ 
tion  for  development  of  material  models.  Specific  issues  were  addressed  concerning  the  dynamic 
response  of  polycrystalline  calcium  carbonate.  This  material  was  chosen  for  this  initial  application 
of  the  experimental  method  to  a  phase- transforming  material  because  extensive  static,  high- 
pressure  and  shock  wave  experiments  had  been  carried  out  previously,9'12  and  the  available  sam¬ 
ples  simplified  execution  and  analysis  of  the  experiments. 

The  experimental  methods  developed  earlier4,5  for  measuring  in-material  particle  velocity  his¬ 
tories  under  combined  compression  and  shear  loading  were  used  with  only  a  few  modifications.  A 
new  experimental  development  was  the  implementation  of  a  suggested  technique13  for  simultane¬ 
ously  measuring  the  two  individual  particle  velocity  components  in  a  single  experiment.  This 
development  permits  precise  time-correlation  between  the  shear  wave  and  the  complex-structured 
compression  wave  that  develops  in  a  phase-transforming  material. 

Time-correlated  shear  and  compression  waves  were  measured  in  the  marble  at  seven  peak 
stress  levels  between  0.5  GPa  and  5  GPa.  The  histories,  measured  at  multiple  depths  within  each 
sample,  included  the  complete  interval  of  compression  and  unloading.  From  these  the  shear, 
compression,  and  longitudinal  release  wave  speeds  were  determined. 

The  shear  waves  were  highly  dispersive  and  were  quenched  by  longitudinal  unloading.  The 
marble  supported  finite  amplitude  shear  waves  over  the  entire  range  investigated  (to  11%  density 
compression  and  4.7  GPa  longitudinal  stress),  which  indicated  that  it  retained  significant  strength. 
Above  3  GPa  longitudinal  stress,  the  shear  wave  amplitude  decreased  substantially  with  propaga¬ 
tion  distance. 

The  longitudinal  waves  had  the  complex  structure  observed  in  previous  measurements,12  14 
and  were  not  self-similar.  The  compression  wave  first  rises  steeply,  and  then  slowly,  forming  a 
knee.  The  continuous,  gradual  rise  after  the  knee  leads  up  to  a  second  steep  portion  of  the  wave. 
The  longitudinal  release  wave  begins  with  a  dispersive  portion  that  leads  down  to  a  sharp  rarefac¬ 
tion  shock  wave.  The  present  measurements  revealed  a  break  in  the  rise  of  the  compression  wave 
for  peak  amplitudes  below  the  knee.  This  break,  and  the  knee  in  the  higher-amplitude 
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experiments,  attenuated  with  propagation  distance.  Hence,  the  compression  waves  were  not  self¬ 
similar. 

A  surface  fitting  technique  was  developed  for  calculating  longitudinal  stress  and  density  in 
unsteady  compression  waves  having  a  wide  variety  of  shapes.  These  calculations  indicated  that 
lack  of  self-similarity  in  the  wave  had  substantial  effects  on  the  stress  and  density  determination 
below  1.55  GPa  longitudinal  stress.  The  locus  of  peak  longitudinal  stress-density  states  was  found 
to  approach  the  hydroetat  for  single  crystal  calcite  for  stresses  above  2.5  GPa,  as  might  occur  if  the 
same  phases  were  produced  under  shock  and  hydrostatic  loading. 

The  shear  wave  speeds  in  shocked  Carrara  marble  as  a  function  of  density  agreed  well  with 
the  ultrasonic  wave  speed  in  hydrostatically  compressed  Oak  Hall  limestone.  The  decrease  in  the 
shear  wave  speed  observed  at  small  compressions  substantiates  earlier  interpretations  based  on 
longitudinal  wave  measurements14’15  that  the  CaCOj(II)  phase  is  produced  under  impact  loading. 
At  higher  compression,  the  shear  modulus  in  the  shocked  marble  increases  with  compression.  This 
occurs  after  the  marble  has  undergone  a  phase  transformation.  This  finding  contradicts  the  infer¬ 
ences  in  earlier  work.14 

The  longitudinal  wave  speeds  in  shocked  Carrara  marble  as  a  function  of  density  were  in  good 
agreement  with  the  speeds  determined  in  Oak  Hall  limestone  under  shock12  and  hydrostatic  load¬ 
ing9  to  9%  density  compression.  At  11%  compression,  the  longitudinal  wave  speed  in  the  shocked 
marble  was  below  an  extrapolation  from  the  hydrostatic  data. 

The  shear  modulus  and  bulk  modulus  in  the  shocked  marble  were  determined  from  the  meas¬ 
ured  wave  speeds  and  calculated  peak  densities.  The  bulk  modulus  in  the  shocked  marble  as  a 
function  of  density  was  in  reasonable  agreement  with  hydrostatic  data  up  to  9%  compression, 
which  is  near  the  onset  on  the  CaC03(III)  phase  stability  field  under  hydrostatic  compression. 
Beyond  this,  at  11%  compression,  there  is  a  pronounced  discrepancy  between  the  bulk  modulus  in 
the  shocked  marble  and  the  bulk  modulus  determined  from  the  hydrostatic  compression  curve  for 
single  crystal  calcite.  No  such  discrepancy  was  indicated  by  the  trend  of  the  peak  stress-density 
states  towards  the  hydroetat.  However,  the  comparison  of  bulk  modulus  data  from  shock  and 
static  experiments  is  direct  and  unambiguous.  The  occurrence  of  a  discrepancy  in  this  comparison 
with  no  corresponding  discrepancy  in  the  comparison  of  the  peak  stress  states  with  the  hydrostat 
indicates  that  the  latter  comparison  must  be  interpreted  with  considerable  caution. 

In  contrast  to  previous  studies,16,17  the  mean  stress-density  relation  for  the  shocked  ma.ble  is 
not  the  integral  of  the  bulk  modulus-density  relation.  Failure  to  account  for  the  systematic  differ¬ 
ence  between  the  frozen  rhase-compoeition  modulus,  which  was  obtained  from  the  wave  speeds  in 
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the  shocked  marble,  and  the  relaxed  modulus,  which  gives  the  slope  of  the  mean  stress-density 
relation,  lead  to  an  overestimate  of  the  mean  stress. 

A  model  of  the  material’s  mean  stress  response  was  presented  to  illustrate  a  method  for  infer¬ 
ring  the  mean  stress-density  relation  for  shocked  materials  that  undergo  phase  transformations 
from  the  experimentally  determined  bulk  modulus-density  relation.  The  model  was  a  variation  of 
the  one  developed  by  Grady  et  al.12  to  describe  the  response  up  through  the  CaC03(II)  phase,  tak¬ 
ing  the  heterogeneity  of  the  stress  field  in  the  polycrystal  into  account. 

Based  on  the  calculated  density,  the  majority  of  the  transformation  occurs  in  the  second  steep 
rise  in  the  compression  wave,  above  a  particle  velocity  of  115  m/s.  The  estimates  of  the  longitudi¬ 
nal  stress  and  density  during  unloading,  suggested  that  the  rarefaction  shock  is  produced  by  direct 
reversion  to  calcite  of  either  CaC03(II)  or  the  succeeding  high-pressure  phase,  depending  on  the 
peak  state  of  compression. 

Comparing  measured  shear  wave  amplitudes  with  previous  wave  propagation  simulations18 
suggested  that  the  marble  exhibits  pressure-dependent  yielding.  This  comparison  indicated  the 
usefulness  of  analyzing  the  shear  wave  amplitudes,  as  well  as  its  complexity.  The  yield  models 
used  in  the  simulations  were  insufficient  to  quantify  the  marble  strength  or  identify  rate  depen¬ 
dence  in  the  yielding. 

Based  on  a  strong  similarity  between  measured  and  simulated18  shear  waves,  it  was  argued 
that  the  compression  wave,  alone,  causes  yielding  in  the  marble  under  combined  compression-shear 
loading  at  a  longitudinal  stress  below  2  GPa.  The  break  in  the  compression  wave  that  occurred 
below  70  m/s  longitudinal  particle  velocity  suggests  that  the  compression  wave,  alone,  may  cause 
yielding  at  only  0.7  GPa  longitudinal  stress. 

The  present  work  has  demonstrated  that  shear  waves  can  be  measured  in  a  shock-induced, 
high-pressure  phase,  providing  additional  information  not  obtainable  from  uniaxial  strain,  plate 
impact  experiments.  The  data  determined  from  analyses  of  the  measured  waves  facilitated  new 
and  more  direct  comparisons  with  static  measurements.  The  difference  in  bulk  modulus  values 
beyond  9%  compression  cannot  be  explained  by  a  difference  between  the  dynamic  mean  stress  at 
which  CaC03(III)  is  produced  and  the  equilibrium  transformation  pressure.  It  suggests  the  onset 
of  a  transformation  to  a  phase  having  a  much  lower  bulk  modulus  than  does  CaCO^(III).  Com¬ 
parison  of  longitudinal  stress-density  states  with  a  hydroetat  is  difficult  to  interpret,  and  was  found 
to  be  misleading  when  judged  against  the  direct  comparison  of  bulk  modulus  data. 
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C.  Time-resolved  optical  spectroscopy  in  shocked  solid s 

The  work  on  fused  silica  and  polycrystalline  calcite  has  clearly  shown  the  value  of  shear  wave 
measurements  to  characterize  the  shocked  state.  However,  there  exists  a  need  for  developing 
experimental  methods  to  probe  the  shocked  state  at  the  microscopic  level.  Hence,  an  important 
objective  of  the  present  work  was  to  explore  the  use  of  optical  spectroscopic  methods  to  obtain 
information  at  the  atomic/molecular  level.  Because  this  is  a  new  undertaking,  the  research  effort 
was  focused  on  developing  experimental  methods  and  performing  feasibility  experiments;  we  have 
emphasized  time-resolved  measurements  in  our  work.  The  equipment  procured  under  AFOSR 
funding  has  been  necessary  to  undertake  this  activity. 

In  one  set  of  experiments  we  developed  an  experimental  method  to  obtain  luminescence  spec¬ 
tra  during  tension.  Ruby  samples  were  shocked  along  the  crystal  c-axis  to  create  successive 
compression  and  tension  loading  using  wave  interactions.  Time-resolved  luminescence  spectra  were 
then  obtained  that  corresponded  to  peak  compression,  complete  unloading  and  peak  tension  states. 
In  tension,  the  luminescence  R  lines  show  a  blue  shift  and  the  line  splitting  increases  with  increas¬ 
ing  tension.  These  data  show  that  the  sapphire  lattice  can  sustain  tensile  stresses  as  high  as  108 
kbar  (under  uniaxial  strain)  without  any  damage.  In  addition,  these  data  have  been  valuable  to 
our  theoretical  effort  to  correlate  optical  spectra  to  symmetry  changes  at  the  local  level.  A 
comprehensive  account  of  this  theoretical  effort  will  be  published  shortly. 

In  a  second  undertaking,  we  developed  a  method  to  do  time-resolved,  Raman  measurements 
in  shocked  solids.  Diamond  was  selected  because  of  its  large  Raman  cross-section19  and  the  consid¬ 
erable  body  of  static  high  pressure  measurements ,20~22  The  frequency  shift  of  the  first  order 
Raman  line  at  uj0  =  1333  cm-1  was  observed  under  uniaxial  strain  along  the  [1 10]  direction  to  a 
longitudinal  stress  of  121  kbar.  The  measured  frequency  shift  was  consistent  with  the  calculated 
strain-induced  changes  in  the  phonon  frequency  in  diamond.  Further  improvements  in  spectral 
resolution  are  needed  to  observe  the  removal  of  the  threefold  degeneracy  for  compression  along  the 
[110]  direction. 

The  present  work  on  time- resolved,  optical  spectroscopy  is  a  good  start  and  makes  a  strong 
case  for  using  this  approach  to  better  understand  shock  deformation  in  solids  in  future  studies. 

D.  Phase  changes  in  calcite  single  crystals  shocked  along  the  x-axis 

This  work  was  motivated  by  the  results  obtained  for  poly  crystalline  CaC03  summarized  ear¬ 
lier.  The  objective  of  the  present  investigation  was  to  undertake  an  exploratory  effort,  using  single 
crystals,  to  examine  the  phase  changes  in  calcite.  Single  crystal  experiments  avoid  the  difficulties 
due  to  stress  concentrations  arising  from  presence  of  randomly  oriented  grains  and  are  necessary  to 
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investigate  detailed  mechanisms.  There  have  not  been  detailed  studies  of  shocked  calcite  crystals 
over  the  compression  range  that  are  comparable  to  measurements  under  hydrostatic  loading.10,11 
Although  it  would  be  desirable  to  use  particle  velocity  gauges,  similar  to  the  polycrystalline  effort, 
available  crystals  did  not  allow  for  such  measurements. 

Quartz  gauges23,24  were  used  to  examine  the  shock  response  of  calcite  crystals  compressed 
along  the  z-axis.  Peak  stress  values  ranged  between  7  and  40  kbar.  Stress  time  profiles  were  meas¬ 
ured  both  at  the  impact  surface  and  after  transmission  through  1-2  mm  thick  samples. 

The  measured  longitudinal  stress-particle  velocity  (<rx  -  u)  response  for  phase  I  is  in  good 
agreement  with  the  calculated  nonlinear  elastic  response  using  elastic  constants  cited  in  the  litera¬ 
ture.11  These  calculations  also  show  the  strong  compressional  anisotropy  in  calcite.  The  mean 
stress-density  curve  for  uniaxial  strain  loading  is  below  the  hydrostatic  loading  curve.  The  calcite 
I — ^11  transition  is  very  rapid  (few  ns  or  less)  and  occurs  at  a  longitudinal  stress  of  approximately 
17.5  ±  0.8  kbar.  The  density  compression  at  the  onset  of  transformation  is  2.1  percent  and  is 
independent  of  the  type  of  loading. 

At  the  higher  stresses  (between  20-40  kbar),  the  shock  response  of  calcite  is  complex  and  the 
transmitted  profiles  show  a  three  wave  structure  with  considerable  rate  dependence  associated  with 
the  second  transformation.  The  stress  at  which  the  second  transformation  occurs  is  considerably 
higher  than  the  pressure  for  the  calcite  II — ^III  transition  under  hydrostatic  loading;  the  reasons  for 
this  difference  are  not  understood.  Because  of  the  time-dependent  features  and  the  impedance 
differences  between  calcite  and  quartz,  a  detailed  analysis  of  the  transmitted  profiles  is  not  feasible. 
The  impact  surface  measurements,  in  which  the  a  —  u  data  are  measured  at  the  calcite/quartz 
interface,  also  reveal  time-dependent  response  beyond  the  calcite  I — ^11  transition.  It  is  not  possible 
to  reconcile  the  impact  surface  and  transmitted  wave  profile  data  with  the  stability  fields  for 
phases  II  and  III  under  hydrostatic  loading.  Our  data  make  a  strong  case  for  the  occurrence  of  a 
final  state  that  is  different  from  the  calcite  III  phase  observed  under  hydrostatic  loading.  These 
findings  support  Aidun’s  results  for  the  polycrystalline  CaCC^. 

Although  our  work  has  not  been  definitive  over  the  22-40  kbar  range,  it  has  provided  interest¬ 
ing  details  and  raised  questions  about  present  understanding  regarding  shock  response  of  calcite. 
The  present  work  suggests  the  need  for  alternate  experiments  for  determining  the  nature  of  the 
shocked  state.  Because  of  the  time-dependent  response,  comparison  of  the  continuum  data  with 
static  measurements  is  difficult  and  the  need  for  spectroscopic  results  to  get  insight  into  the  crystal 
structure  under  shock  loading  is  emphasized. 
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m,  SUMMARY  OF  THE  LASER  ENERGY  DEPOSITION  EFFORT 


The  investigation  of  non-linear  laser  pulse  interaction  with  transparent  solids  culminated  in 
explaining  the  fundamental  mechanisms  responsible  for  energy  deposition  and  the  material's 
response  to  it.  As  a  result  we  have  provided  an  explanation  for  single-shot  and  multi-pulse  intrin¬ 
sic  laser  damage  in  wide-gap  high  laser  power  optical  materials. 

In  the  course  of  this  work  new  experimental  techniques  have  been  developed  for  the  measure¬ 
ment  of  energy  deposition  from  a  laser  pulse  at  visible  and  near  infrared  wavelengths  into  tran¬ 
sparent  solids.  Their  application  to  alkali  halides  and  fused  silica  (Si02)  revealed  that  lattice  heat¬ 
ing  occurs  via  nonlinear  absorption  of  laser  photons  by  multiphoton  excited  free  electrons  with  only 
small  contributions  that  can  be  attributed  to  point  defects  created  during  the  laser-solid  interac¬ 
tion. 

In  contrast  to  these  findings  is  our  discovery  that  multi-pulse  effects  at  intensities  below  the 
single-shot  damage  threshold  are  dominated  by  laser-generated  stable  defects.  We  have  formulated 
a  new  model  of  intrinsic  multi-pulse  bulk  damage  that  is  consistent  with  our  experimental  observa¬ 
tions.  It  is  based  on  thermomechanical  stress  induced  by  accumulation  of  multiphoton  generated 
F-centers. 

The  interaction  of  intense  laser  pulses  with  optical  materials  has  been  the  subject  of  study  for 
over  two  decades  with  the  goal  to  explain  the  basic  mechanisms  of  laser  damage  (optical  dielectric 
breakdown)  of  intrinsically  pure  and  transparent  wide-gap  solids.  Alkali  halides  have  traditionally 
been  of  particular  interest  in  these  investigations. 

In  discussing  intrinsic  optical  dielectric  breakdown  in  alkali  halides  one  has  to  distinguish 
between  single-pulse  damage  and  catastrophic  failure  due  to  repetitive  exposure  of  the  same 
material  to  laser  pulses  of  power  density  below  the  single  pulse  damage  threshold.  In  both  cases 
breakdown  is  the  result  of  energy  deposition  from  the  laser  beam  to  the  solid  associated  with  a 
sequence  of  complex  electronic  and  atomic  (photo-chemical)  processes,  but,  as  we  will  show  in  this 
report,  with  significantly  different  importance  of  laser-induced  defects. 

The  primary  process  in  both  cases  is  electron-hole  pairs  (exciton)  formation  by  multiphoton 
interband  transitions.1  Single  pulse  damage  is  due  to  non-linear  absorption  by  the  resulting  free 
electrons  with  only  small  contributions  from  absorption  by  primary  point  defects  that  are  produced 
as  follow-on  products  of  exciton  recombination.  In  contrast  to  this  situation,  the  shot  to  shot  accu¬ 
mulation  of  stable  point  defects  and  the  associated  thermo-mechanical  stresses  turned  out  to  be  the 
cause  of  multi-pulse  failure.1 

Historically,  the  mechanism  of  intrinsic  laser  damage  was  thought  to  be  avalanche  impact  ion¬ 
ization  by  free  electrons  which  had  gained  sufficient  kinetic  energy  during  the  interaction  with  the 
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intense  laser  photon  field.2  The  origin  of  this  concept  can  be  traced  to  the  work  by  von  Hippel3 
and  others  who  believed  dc  dielectric  breakdown  occurs  by  this  process.  Forty  years  later,  Yablo- 
novitch  and  Bloembergen  noted  an  apparent  agreement  of  laser  damage  field  strengths  (all  smaller 
than  20  MV/cm)  at  10.6  /im  and  dc  data.  This  led  them  to  propose  an  early  version  of  the  so- 
called  electron  avalanche  laser  breakdown  theory4  which  until  very  recently  dominated  the  think¬ 
ing  of  virtually  every  worker  in  the  field.5 

The  picture  of  dielectric  breakdown  at  dc  to  optical  frequencies  has  dramatically  changed  in 
the  last  five  years  due  to  the  work  of  Fischetti,  DiMaria,  Cartier,  and  co- workers  at  IBM  and  our 
group  at  Washington  State  University.  Cartier  et  al.6  were  able  to  measure  the  electron-photon 
scattering  rates  in  amorphous  SiC>2  as  a  function  of  electron  energy,  confirming  calculations  by 
Fischetti  et  al.  and  enabling  Monte  Carlo  simulations  at  high  dc  field  of  electron  transport  in  a 
number  of  wide-gap  solids.7  DiMaria  and  co-workers  developed  methods  for  the  direct  measure¬ 
ment  of  electron  energy  distributions  under  high  field  conditions.8  In  agreement  with  Fischetti’s 
predictions  they  found  that  the  energy  distributions  provided  no  evidence  for  electron  multiplica - 
tion  (avalanche  formation)  up  to  dc  field  strengths  of  20  MV/cm.9  It  turned  out  that  efficient 
scattering  of  electrons  above  2  eV  with  acoustic  phonons  stabilizes  the  electron  energy  distribution 
and  prevents  runaway.10  Similar  conclusions  are  expected  to  hold  for  alkali  halides. 

Independent  of  the  fundamental  investigations  at  IBM  we  reached  similar  conclusions  con¬ 
cerning  intrinsic  damage  in  NaCl,  KBr,  KI,  and  Si02  at  optical  frequencies.  We  were  able  to 
develop  experimental  techniques  for  the  measurement  of  prebreakdown  energy  deposition  from  the 
laser  pulse  into  intrinsically  transparent  solids.  These  made  possible  the  systematic  investigation 
of  the  processes  involved  in  the  deposition  of  energy  from  the  laser  beam  into  the  transparent  solid 
(including  the  role  of  laser-generated  primary  defects)  and  provided  hard  evidence  for  the  absence 
of  avalanche  formation  at  visible  and  near-infrared  wavelengths  up  to  the  intrinsic  single-shot  dam¬ 
age  threshold. 

A.  The  mechanism  of  single  pulse  laser  damage  in  alkali  halides 

Two  types  of  experiments  were  performed  in  the  prebreakdown  regime  to  measure  the 
absorption  of  laser  beam  energy,  and  these  measurements  were  compared  with  models  including 
processes  that  are  known  to  occur  in  alkali  halides  following  electron-hole  pair  production.  Two 
alternative  models  for  free-electron  heating  were  investigated.  The  difference  in  (theoretical)  heat¬ 
ing  efficiency  of  these  two  lead  to  conclusions  regarding  the  relative  importance  of  primary  defects 
in  single-shot  laser  damage. 
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Calibrated  photoacoustic  measurements  of  prebreakdown  energy  absorption  from  532  nm 
laser  pulses  in  NaCl  revealed  that  the  primary  process  for  the  interaction  is  four-photon  absorption. 
These  results  are  compared  with  three  models  of  the  interaction. 

They  account  in  detail  for  the  creation  and  recombination  of  e-h  pairs,  creation  and  annihila¬ 
tion  of  primary  defects  (Vk,  STE,  STH,  F),  and  the  resulting  heating  of  the  lattice.  We  allow  for 
two  possible  pair  generation  processes,  multiphoton  absorption  and  impact  ionization.  Generally, 
electrons  are  elevated  to  the  conduction  band  either  from  the  valence  band  or  from  ionized  defects 
(STE’s,  F-centers).  Free  valence  band  holes  are  trapped  rapidly,  but  may  be  freed  again  by  hole 
absorption,  then  quickly  retrapped.  The  net  result  is  heat  deposition  of  ~hj  per  event.  Vk  centers 
may  orient  themselves  into  directions  relatively  inert  to  the  polarized  laser  light,  thus  the  hole 
absorption  does  not  occur  for  each  Vk  center  throughout  the  duration  of  the  laser  pulse. 

F-center  formation  is  assumed  to  occur  through  the  triplet  STE  state  at  high  temperatures. 
This  effect  is  included  or  excluded  as  to  its  interest  in  the  particular  experiment  modeled. 

The  two  treatments  of  free-carrier  heating  are  called  (1)  the  free-electron  model,  a  nonlinear 
absorption  of  electromagnetic  energy  by  fairly  energetic  conduction  electrons  and  (2)  the  polaron 
model,  wherein  low  energy  conduction  electrons  which  are  strongly  coupled  to  lattice  ions  absorb 
photons  linearly.  The  first  is  vastly  more  efficient  theoretically,  but  there  is  no  a  priori  reason  to 
prefer  one  of  these  unproven  models.  For  completeness,  we  also  compared  the  data  with  a  model 
including  electron-impact  ionization  avalanche  production  of  e-h  pairs. 

Both  treatments  of  free  carrier  heating  satisfactorily  represent  the  measurements.  The  four- 
photon  absorption  cross-section  is  ambiguous  because  the  photoacoustic  effect  is  a  result  of  total 
absorbed  energy.  Allowing  for  secondary  absorption  processes  (in  addition  to  four-photon  absorp¬ 
tion)  without  a  preference  for  one  of  the  two  alternatives,  it  cannot  be  determined  which  value  of 
<7^  is  most  correct,  since  both  are  satisfactory.  If  model  (1)  is  correct,  then  fewer  electrons,  and 
therefore,  fewer  primary  defects  are  required  to  produce  the  effect.  In  fact,  for  model  (l),  »  92% 
of  all  energy  absorption  is  via  free  electron  heating.  In  model  (2)  defect  absorption,  primarily  Vk 
center  absorption,  accounts  for  most  of  the  heating  (>  50%).  Further,  the  avalanche  process  does 
not  account  for  the  data. 

In  order  to  discover  which  model  it  most  correct,  and  thus  the  importance  of  primary  defects, 
it  is  necessary  to  model  similar  measurenents  with  a  known  value  of  .  This  was  accomplished 
using  the  self-trapped  exciton  recombination  luminescence  or  STERL  technique. 

The  STERL  technique  of  measuring  multiphoton  excitation  takes  advantage  of  this  intrinsic 
radiative  relaxation  of  electron-hole  (e-h)  pairs  in  alkali  halides.  Since  the  process  is  intrinsic,  the 
total  light  signal  output  is  a  measure  of  the  number  of  e-h  pairs  generated  in  an  excitation  event. 
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This  then  represents  an  improvement  in  measuring  multiphoton  cross-sections  over  the  photoacous¬ 
tic  method  because,  in  a  limited  range  of  incident  pulse  intensities,  the  response  of  the  sample  to 
incident  ionizing  radiation  is  proportional  only  to  the  total  number  of  excited  carriers.  This  pro¬ 
portionality  falls  off  when  the  crystal  temperature  exceeds  that  at  which  non-radiative  recombina¬ 
tion  processes  become  important.  This  effect  was  used  as  a  thermometer  to  determine  lattice  heat¬ 
ing  in  KBr  when  pulse  energies  were  increased  beyond  the  level  where  thermal  quenching  of  the 
luminescence  begins.  Armed  with  the  value  of  obtained  at  low  pulse  intensities,  the 
multiphoton-polaron  and  -free  electron  theories  of  energy  absorption  were  compared  with  experi¬ 
mented  data,  revealing  that  the  polaron-defect  absorption  is  simply  too  weak  to  account  for  the 
observed  temperature  increase,  leaving  the  multiphoton-free-carrier  model  as  a  complete,  satisfac¬ 
tory  model  of  prebreakdown  energy  absorption.  Further,  temperature  data  indicate  that  the 
highest  intensity  nondamaging  laser  pulses  induced  heating  in  KBr  to  only  a  small  amount  below 
the  melting  point.  These  measurements  could  all  be  explained  and  modeled  without  any  necessity 
of  including  avalanche  generation.  In  fact,  these  measurements  and  calculations  preclude  the  role 
of  avalanche  generation  at  least  up  to  the  melting  point,  even  though  electron  densities  greater 
than  1018cm  3  are  present.  These  results  were  verified  in  KI  (three  photon  band  gap  at  532  nm). 

The  conclusion  reached  in  the  STERL  study  was  bolstered  by  a  third  experiment.  Pulses  of 
266  nm  laser  light  were  used  to  create  e-h  pairs  in  NaCl  and  Si02  (via  2-photon  absorption),  and 
the  photoacoustic  signal  was  monitored  as  a  function  of  coincident  pulse  energies  of  1064  nm  laser 
light.  Since  this  wavelength  is  too  long  to  create  e-h  pairs,  the  only  mechanism  by  which  absorp¬ 
tion  could  take  place  is  via  conduction  carriers.  The  nonlinear  nature  of  the  resultant  signal  veri¬ 
fied  the  correctness  of  the  free-electron  treatment  of  the  conduction  carrier  absorption,  thus 
relegating  the  participation  of  defects  to  a  minor  or  possible  insignificant  role  in  the  single  hot  dam¬ 
age  process.  Exaggerating  their  roles  in  the  model  calculations  showed  that  the  models  would  cease 
to  represent  the  data.  Therefore,  we  conclude  primary  defects  play  at  most  a  minor  role  in  single 
pulse  laser  damage  of  pure  alkali  halides. 

The  work  presented  up  to  this  point  reflects  significant  advances  in  the  understanding  of  the 
fundamental  interactions  between  high  power  pulsed  laser  beams  and  transparent  solids.  Identify¬ 
ing  the  intrinsic  single  shot  damage  threshold  provides  an  important  standard  of  performance  for 
optical  materials.  But  perhaps  of  greater  importance  from  a  practical  standpoint  are  the  processes 
contributing  to  multi-pulse  laser  damage.  In  the  next  section,  we  explore  a  new  model  of 
multiple-pulse-on-one-site  damage  based  on  the  accumulation  of  laser  induced  color  center  forma¬ 
tion.  This  new  model,  an  out-growth  of  the  model  calculations  discussed  in  previous  sections,  is 
not  concerned  with  electrons  but  in  the  alteration  of  the  mechanical  properties  of  the  crystalline 
solid  which  accompany  the  formation  of  lattice  defects. 


-14- 


B.  Multiple  shot  bulk  damage 

In  contrast  to  single  shot  experiments,  where  temperature  increases  of  several  hundred 
degrees  were  observed,1  during  the  multiple  pulse  experiments  presented  here  the  peak  temperar 
ture  increase  is  at  most  a  few  tens  of  degrees.  Under  these  conditions  it  is  clear  that  a  different 
failure  mode  is  required  to  explain  the  eventual  damage  of  a  crystal  at  laser  power  levels  well  below 
the  single  shot  damage  threshold.  The  most  obvious  cumulative  effect  which  might  suffice  to  cause 
damage  at  power  levels  below  the  single  shot  damage  threshold  is  the  buildup  of  F-centers  as  a 
result  of  laser- induced  electron-hole  pair  generation. 

In  alkali  halides,  F-centers  have  a  typical  absorption  peak  of  2-3  eV  and  can,  therefore, 
release  their  electrons  into  the  conduction  band  through  single  or  two-photon  absorption  at  visible 
frequencies.  However,  Jones  et  al.1  have  shown  that  in  alkali  halides  the  effect  of  F-centers  on  lat¬ 
tice  heating  is  insignificant,  even  with  relatively  large  initial  concentrations.  Thus,  at  first,  it  would 
appear  that  F-center  accumulation  is  not  responsible  for  the  observed  multi-shot  damage  behavior 
in  KBr  and  KI  reported  here.  However,  we  now  suggest  an  alternative  mechanism  in  which  these 
defects  can  make  major  contribution  in  the  damage  process. 

It  has  been  shown  in  alkali  halides  that  when  a  crystal  is  radiatively  colored  it  undergoes  an 
expansion  due  to  the  production  of  Frenkel  defects.  As  F-centers  accumulate  in  the  focal  volume 
from  shot  to  shot,  local  stresses  will  develop  since  the  associated  expansion  is  impeded  by  the  sur¬ 
rounding  uncolored  crystal.  If  this  stress  reaches  a  sufficient  level  it  is  expected  that  the  material 
will  crack  or  experience  plastic  flow  resulting  in  a  nonuniformity  in  the  crystal  which  acts  as  a 
scattering  center.  Scattering  of  the  incident  laser  pulse  from  this  nonuniformity  is  detected  and 
interpreted  as  damage. 

The  model  presented  here  assumes  that  failure  occurs  when  the  local  stress,  a,  exceeds  the 
yield  stress,  cry.  To  simplify  matters  we  take  the  local  stress  to  be  the  mean  of  the  diagonal  ele¬ 
ments  of  the  stress  tensor,  or  a  =(<xu  +  a?!  +033) /3.  The  yield  stress,  however,  is  not  indepen¬ 
dent  of  defect  concentration.  Alkali  halide  crystals  exhibit  an  increase  in  the  yield  stress  after  radi¬ 
ative  coloring.  Nadeau  et  al.11  found  the  increase  in  yield  stress  to  be  proportional  to  the  square 
root  of  the  F-center  concentration,  np,  or 

<?y  =  A  +  By/n p" ,  (1) 

with  A  =  6.47  x  106  dyne  cm"2  and  B  =  6.54  x  10"2  dyne  cm-1/2  for  KBr. 

The  induced  stress  consists  of  two  components  of  different  origin.  The  first  is  due  to  the  lat¬ 
tice  expansion  in  the  focal  volume  caused  by  the  production  of  F-centers.  The  associated  volume 
increase,  u,  due  to  the  production  of  one  F-center  has  been  determined  by  measuring  the  linear 
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expansion  of  a  crystal  as  a  function  of  F-center  concentration:  v  =  0.074a3  for  NaCl  and  KBr  (a3 
is  the  volume  of  the  original,  undisturbed  lattice  cell). 

In  calculating  the  thermal  stress,  <7x,  we  are  able  to  make  use  of  our  single  shot  results  which 
showed  that  the  peak  temperature  increase,  during  the  laser  pulse,  is  very  nearly  proportional  to 
Fq1,  where  F0  is  the  peak  photon  flux  density,  and  m  is  the  order  of  multiphoton  absorption. 
Therefore,  the  thermal  stress  is 

*  CaKFg1.  (2) 

Here,  a  is  the  coefficient  of  thermal  expansion,  C  is  a  constant  of  proportionality  and  F0  is  the 
peak  photon  flux.  The  total  stress  is  the  sum  of  these  two  components 

o  =  op  4*  0“t*  (3) 


It  remains  to  determine  the  growth  of  the  F-center  concentration.  Drawing  an  analogy  to  the 
failure  of  materials  under  load  Zhurkov  et  al.12  and  Manenkov  et  al.13  have  suggested  that  multi¬ 
shot  damage  may  be  due  to  a  mechanochemical  reaction  governed  by  a  reaction  rate  of  the  form, 
dN/dt  «  kbT  exp{  -(U0  -  7^) A  bT}.  N  represents  the  population  of  the  appropriate  reaction  pro¬ 
duct  associated  with  the  damage  process,  kB  is  Boltzmann’s  constant,  T  is  the  temperature,  and  U0 
is  the  activation  energy  of  the  process,  a  is  then  the  applied  stress  and  7  is  a  material  dependent 
parameter.  The  effect  of  the  stress  is  to  reduce  the  activation  energy  for  the  nucleation  process. 
We  propose  that  F-center  production  is  the  pertinent  mechanochemical  process  to  which  the  above 
rate  applies.  Indeed,  there  have  been  studies  where  an  increase  in  F-center  production  with 
applied  stress  was  observed. 

The  number  of  F-centers  produced  during  any  one  laser  pulse  under  stress  free  conditions  is  a 
known  function  of  temperature.  F-center  production  is  also  proportional  to  the  number  of 
electron-hole  pairs  produced.  For  KBr  and  KI  at  532  nm,  this  is  in  turn  proportioned  to  F™. 
Thus,  for  any  one  laser  pulse 


AnF  =  DF“kBTexp<  -  U0  -  7<r  A  bT 


(4) 


with  D  being  a  constant  of  proportionality,  which  is  determined  from  the  known,  stress  free,  defect 
production  efficiencies.  An  iterative  calculation  is  then  used  to  follow  the  buildup  of  F-centers  and 
stress  from  one  shot  to  the  next  until  the  condition  for  failure  is  met. 

The  experimental  procedure  is  to  expose  the  sample  to  a  series  of  laser  pulses  of  uniform  peak 
power  at  a  repetition  frequency  of  about  one  Hertz.  The  laser  pulses  are  Gaussian  in  both  spatial 
and  temporal  profiles,  with  a  spot  size  of  8-10  jim  and  a  pulse  length  of  approximately  70  psec. 
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These  values  refer  to  the  diameter  and  half-width  at  1/e  of  the  peak  photon  flux  density,  respec¬ 
tively.  The  peak  flux  density  varied  from  about  0.9  x  1029  to  1.6  x  10 29  photons  cm-2  sec-1  (34-52 
GW/cm2)  for  the  experiments  on  KBr  and  ranged  from  0.6  x  1029  to  1.0  x  1029  photons  cm”2 
sec-1  (22-37  GM/cm2)  for  KI.  The  sample  is  mounted  on  the  cold  finger  of  a  closed  cycle,  optical 
cryostat  and  the  temperature  can  be  varied  from  about  50  K  to  room  temperature.  The  condition 
of  the  crystal  is  monitored  by  measuring  the  scattering  of  laser  light  perpendicular  to  the  beam  axis 
with  a  photomultiplier  tube.  Damage  is  assumed  to  have  occurred  when  a  sudden  increase  in  light 
scattering  is  observed.  The  decay  time  for  the  thermal  disturbance  is  on  the  order  of  10”5  sec  for 
the  focal  spot  size  used  in  the  experiment,  thus  precluding  a  simple  buildup  of  temperature  at  the 
pulse  repetition  frequency  used. 

It  is  important  to  differentiate  between  data  which  represent  the  intrinsic  behavior  of  the 
material  and  those  which  are  extrinsic.  The  intrinsic  behavior  represents  a  fundamental  limit  of 
the  crystal,  whereas  the  extrinsic  behavior  can  be  caused  by  any  number  of  defects  or  impurities 
initially  present  in  the  crystal.  Although  for  any  give  photon  flux  density  there  may  be  a  large 
variation  in  the  number  of  shots  withstood  by  different  sites  in  the  sample,  all  data  points  are  seen 
to  lie  within  a  distinct  envelope.  Rather  than  performing  a  statistical  analysis  of  the  variations  we 
feel  that,  given  a  sufficient  number  of  points,  the  intrinsic  behavior  of  the  crystal  becomes  evident. 
The  intrinsic  limit  exhibits  itself  as  a  boundary  beyond  which  no  data  points  fall. 

This  procedure,  however,  can  only  be  applied  to  crystals  of  very  high  purity.  Otherwise, 
there  may  be  too  few  sites  which  exhibit  intrinsic  behavior  for  the  envelope  of  points  to  become 
clear.  For  this  reason  we  use  crystals  from  the  University  of  Utah  grown  with  reactive  atmosphere 
processes,  ultra  pure  material. 

At  low  initial  temperatures,  where  the  efficiency  of  F-center  production  is  very  small,  the 
thermal  stress  is  particularly  important.  In  this  case  the  stress  must  reach  a  sufficient  level  to 
boost  the  F-center  formation  efficiency  to  start  the  accumulation  process,  and  we  see  a  distinct 
threshold-like  behavior.  At  higher  temperatures  the  F-center  formation  rate  is  large  enough 
(without  being  enhanced  by  the  thermal  stress)  so  that  there  is  a  much  weaker  flux  dependence. 

Results  of  single  shot  experiments  indicate  that  defects  play,  at  most,  a  minor  role  in  the 
heating  (and  thus  damage)  of  wide  bandgap,  optical  materials.  However,  at  the  power  levels  used 
in  multishot  experiments,  heating  alone  is  not  sufficient  to  cause  damage.  In  this  case  the  accumu¬ 
lation  of  defects  form  shot-to-shot  leads  to  a  buildup  of  stress  and,  ultimately,  the  failure  of  the 
material.  We  have  developed  a  model  for  the  behavior  of  KBr  and  KI  exposed  to  multiple  laser 
pulses  induced  in  the  crystal.  The  model  applies  to  the  intrinsic  behavior  of  the  crystal  in  that  it 
does  not  appeal  to  a  probabilistic  mechanism  or  to  the  inclusion  of  absorbing  impurities  or  defects. 
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Instead,  it  attempts  to  incorporate  known  responses  of  the  crystal  to  the  laser  pulses,  temperature 
and  stress.  At  this  stage  in  its  development  this  model  predicts  the  correct  dependence  of  multi¬ 
shot  laser  damage  on  temperature  and  photon  flux  density  for  the  two  materials  studied.  However, 
additional  refinements  to  the  model  (such  as;  the  inclusion  of  strain  rate  effects  and  the  possibility 
of  defect  clustering)  and  experimental  confirmation  of  the  defect  accumulation  as  described  in  the 
model  are  required  before  a  final  assessment  of  its  validity  can  be  made. 
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IV.  CONCLUDING  REMARKS 


The  shock  wave  effort  has  demonstrated  the  usefulness  of  making  shear  wave  measurements 
for  understanding  the  nonlinear  response  of  materials  under  dynamic  loading.  In  fused  silica,  the 
shear  wave  measurements  led  to  the  finding  of  reversible,  shea’*  enhanced  companion.  In  addition, 
stress  de viators  determined  from  our  work  are  considerably  different  from  those  assumed  in  present 
studies. 

The  work  on  polycrystalline  calcite  represents  the  first  study  to  use  shear  wave  measurements 
to  understand  shock  induced  phase  changes.  The  shear  wave  measurements  have  provided  new 
information  about  the  calcite  response  that  contradicts  existing  infer  nces  about  the  shock  response 
in  calcite.  The  direct  measurement  of  shear  modulus  decrease,  associated  with  the  calcite  I— >11 
transition,  is  an  important  result  that  confirms  the  occurrence  of  this  transition.  More  importantly, 
the  subsequent  increase  in  shear  modulus  at  higher  compressions  was  a  surprising  result.  This 
measurement  underscores  the  need  to  quantify  the  shocked  state  in  more  detail  and  to  use  caution 
in  comparing  the  ax  -  V  results  under  shock  loading  with  static  P  -  V  data  to  understand  the 
shock  response.  The  possibility  that  a  phase  other  than  calcite  III  is  being  formed  under  shock 
loading  is  an  interesting  one  and  that  needs  to  be  pursued  in  future  studies.  More  definitive  con¬ 
clusions  will  require  optical  spectroscopy  measurements  in  conjunction  with  use  of  calcite  single 
crystals.  The  small  effort  undertaken  on  the  shock  measurements  in  calcite  single  crystals  reveals 
clearly  the  I— >11  transition.  However,  the  response  between  22-40  kbar  is  complex  and  character¬ 
ized  by  time  dependent  transformation  to  at  least  one  or  more  phases.  Continuum  measurements 
alone  are  insufficient  to  resolve  this  response. 

During  the  course  of  the  present  work,  we  demonstrated  the  ability  to  do  time-resolved,  opti¬ 
cal  spectroscopy  experiments  in  shocked  solids.  Both  luminescence  and  Raman  measurements  can 
now  be  carried  out.  These  developments  represent  new  approaches  to  characterizing  the  shocked 
state  and  are  complementary  to  the  continuum  measurements.  The  combination  of  the  two  types 
of  measurements  will  be  required  to  better  understand  tne  shocked  state. 

The  laser  effort  was  successful  in  clearing  up  several  confused  issues  in  the  field  of  interaction 
of  intense  laser  pulses  with  intrinsically  transparent  optical  materials.  In  fact,  the  physical 
processes  responsible  for  laser  damage  due  to  energy  deposition  from  the  laser  pulse  into  the  tran¬ 
sparent  solid  were  explained  by  carefully  executed  measurements  in  the  prebreakdown  power  den¬ 
sity  regime.  As  a  result,  we  have  shown  that  the  presumption  of  laser  damage  due  to  laser-induced 
electron  avalanche  formation  is  not  needed  to  explain  intrinsic,  single  pulse  bulk  damage  for 
Eg/&j<5  (Eg  is  the  width  of  the  forbidden  gap  of  the  solid)  and  that  avalanches  are  not  measur¬ 
ably  formed  in  NaCl  and  Si02  for  Eg/ftu;>7  even  when  large  concentrations  of  free  electrons  are 
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present  and  the  rms  electric  field  is  near  the  breakdown  value.  Instead,  the  process  at  visible 
wavelengths  is  satisfactorily  explained  by  laser  energy  absorption  by  multi-photon-excited  free  elec¬ 
trons.  Damage  correlates  with  lattice  heating  to  near  the  melting  point. 

Based  on  our  results,  we  can  define  the  intrinsic,  single-shot,  bulk  damage  for  three  alkali 
halides  at  A  =  532  nm,  r  —  80-100  psec 


KI 

:  =  (6  ±  5)  x  10  8Icm6sec2 

IB  =  75  ±  7.5  GW/cm  2 

KBr 

:  a<4>  =  (2  ±  1)  X  10_u2cm8sec3 

IB  =  240  ±  24  GW/cm  2 

NaCl 

:  <r(4)  =  1.5  x  10~114cm8sec3 

IB  =  745  ±  75  GW/cm  2 

The  reported  errors  are  simple  uncertainties  in  pulse  energy,  beam  radius  and  pulse  length. 

Within  the  multiphoton  initiated  damage  model,  dilute  impurities  or  defects  with  interband 
electronic  states  do  not  make  a  significant  contribution  to  single  pulse  damage  in  terms  of  energy 
absorption  or  providing  free  electrons.  However,  laser-induced  lattice  defect  formation  at  intensi¬ 
ties  below  the  single  shot  damage  threshold  appears  to  be  a  plausible  explanation  of  multipulse  opt¬ 
ical  damage.  This  new  avenue  of  investigation  leads  away  from  the  idea  that  defects  provide 
avalanche  initiating  electrons,  thus,  simply  lowering  the  single  shot  threshold. 

Using  the  knowledge  of  single  pulse  interactions  as  a  starting  point,  we  have  developed  an 
explanation  for  multishot  bulk  damage  based  on  the  response  of  the  crystal  to  the  laser  induced 
temperature  distribution  and  the  accumulation  of  stable  defects.  The  experimental  techniques  we 
have  developed  have  been  met  with  great  interest  by  the  so-called  laser  damage  community  which 
convenes  annually  at  Boulder,  CO  for  the  Laser  Damage  Symposium.  They  were  transferred  to 
the  Advanced  Drivers  Program  of  the  Internal  Confinement  Fusion  Effort  and  the  Y-Division 
(Lasers)  at  the  Lawrence  Livermore  National  Laboratory  (Contacts:  Drs.  Mark  R.  Kozlowski  and 
C.  Robert  Wolfe).  The  intended  applications  there  are  in  the  study  of  laser  pulse  interaction  with 
dielectric  coatings  on  high  power  laser  optical  components. 

Our  investigations  of  the  response  of  transparent  optical  materials  to  rapid  energy  deposition 
from  intense  laser  pulses  created  a  firm  experimental  basis  for  the  understanding  of  the  fundamen¬ 
tal  processes  involved  in  dielectric  breakdown  at  optical  frequencies.  This  understanding  led  away 
from  the  traditional  picture  of  avalanche  breakdown,  and  the  agreement  of  our  work  with  recent 
new  conclusions  concerning  dc  breakdown  is  very  gratifying. 

We  conclude  by  noting  that  understanding  the  dynamics  of  the  laser  damage  interaction  in 
multiple  shot  experiments  will  require  time-resolved  optical  measurements  to  quantify  the  thermo- 
mechanical  stresses  around  the  defects.  The  time- resolved,  optical  spectroscopy  measurements 
under  the  well  defined  conditions  of  plane  shock  wave  loading  are  a  start  in  this  direction. 
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Ruby  Crystals" 

11:15:  Discussion 

12:00:  No  Host  Luncheon  at  the  Compton  Union  Building 

Laser  Damage  o/  Optica/  Materials 

1:30:  P.F.  Braunlich,  "Historical  Review  of  Laser  Damage  Studies:  Theories,  Methods,  Contro¬ 

versies" 

2:00:  S.C.  Jones,  "Models  for  and  Measurements  of  Multiphoton  Initiated  Prebreakdown 

Energy  Absorption  in  Wide  Gap  Insulators' 
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Laser  DamaQc  of  Optical  Materials  —  Continued 


2:30:  X.A.  Shen,  "Multiphoton  Absorption,  Lattice  Heating  and  Laser  Damage  in  Alkali 

Halides" 

3:00:  Break 

3:15:  R.T.  Casper,  "Models  and  Measurements  for  Single-  and  Multi-pulse  Laser  Damage  in 

Alkali  Halides" 

3:45:  Discussion 

Thursday,  May  18 

Invited  Lectures : 

8:15:  J.W.  Shaner,  "Soft  Sphere  Modeling  of  Fluid  Phase  Thermodynamics:  Applications  to 

Rapid  Energy  Deposition" 

9:05:  Break 

9:25;  B.S.  Holmes,  "An  Overview  of  Recent  Laser  Effects  Studies  at  SRI" 

10:15:  Discussion  and  Closing  Remarks 

10:45:  Tour  of  the  Shock  Dynamics  Laboratory  and  the  Laser  Physics  Laboratory. 
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